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We develop an interferometric technique for making time-resolved measurements of field- 
quadrature operators for nonequilibrium ultracold bosons in optical lattices. The technique exploits 
the internal state structure of magnetic atoms to create two subsystems of atoms in different spin 
states and lattice sites. A Feshbach resonance turns off atom-atom interactions in one spin sub- 
system, making it a well-characterized reference state, while atoms in the other subsystem undergo 
nonequilibrium dynamics for a variable hold time. Interfering the subsystems via a second beam- 
splitting operation, time-resolved quadrature measurements on the interacting atoms are obtained 
by detecting relative spin populations. The technique can provide quadrature measurements for 
a variety of Hamiltonians and lattice geometries (e.g., cubic, honeycomb, superlattices), including 
systems with tunneling, spin-orbit couplings using artificial gauge fields, and higher-band effects. 
Analyzing the special case of a deep lattice with negligible tunneling, we obtain the time evolu- 
tion of both quadrature observables and their fluctuations. As a second application, we show that 
the interferometer can be used to measure atom-atom interaction strengths with super-Heisenberg 
scaling n" 3 / 2 in the mean number of atoms per lattice site, and standard quantum limit scaling 
M ~ x l 2 in the number of lattice sites. In our analysis, we require M 3> 1 and for realistic systems 
n is small, and therefore the scaling in total atom number N = nM is below the Heisenberg limit; 
nevertheless, measurements testing the scaling behaviors for interaction-based quantum metrologies 
should be possible in this system. 

PACS numbers: 37.10.Jk, 37.25. +k, 67.85.Hj, 06.20.Dk 



Ultracold atoms in optical lattices [IH5] are versa- 
tile systems for studying nonequilibrium physics [3Mll|. 
particularly for making time-resolved measurements of 
system observables. For example, lattice collapse-and- 
revival experiments [4j IT2HI4] follow atom populations in 
the k = quasi-momentum of nonequilibrium Bose gases 
versus system evolution time. Our goal is to extend the 
power of these experiments to enable time-resolved mea- 
surements of a greater variety of system observables and 
Hamiltonians. A second goal is to use nonequilibrium 
collapse-and-revival dynamics to investigate the physics 
of interaction-based (or nonlinear) quantum metrology, 
which exploits multibody dynamics to characterize a sys- 
tem and its parameters [T5Tl22| . 

In this Letter, we design an interaction-based in- 
terferometer for making time-resolved measurements 
of quadrature operators |23| of nonequilibrium matter 
fields in optical lattices. Measurement of quantum-field 
quadratures X k = (e~ lC - k A k + e li > k A\)/2, where A k anni- 
hilates atoms with quasi-momentum k and spin state a, 
and Cfe is a tunable phase, can provide direct determina- 
tion of the order parameter for a Bose gas [23] , as opposed 
to inferring it from density measurements [25, 26 . The 
technique, analogous to homodyne detection, works as 
follows. Exploiting the internal state structure of ultra- 
cold atoms, we split a supcrfluid into two subsystems-the 
arms of the interferometer-consisting of atoms in differ- 
ent spin states and lattice sites. A Feshbach resonance 



[27] turns off atom-atom interactions in one spin sub- 
system, making it a well-characterized reference "beam" , 
while atoms in the other subsystem undergo nonequilib- 
rium dynamics for a variable hold time. Interfering the 
subsystems via a second beam-splitting operation, spa- 
tial quadrature observables are obtained by measuring 
spin populations. 

As a second, distinct application of the interferome- 
ter, we show that atom-atom interaction strengths can be 
determined from the frequencies of matter-wave collapse- 
and-revival oscillations with "super-Heisenberg" [T5T - fT5] 
scaling fi~ 3 ' 2 in the mean number of atoms per lattice site 
n, and standard quantum limit (SQL) scaling M 1 / 2 in the 
number of lattice sites M. In contrast, we find that the 
optimal scaling from conventional collapse-and-revival 
dynamics, which do not exploit entanglement of motional 
and spin degrees of freedom, is n~ 3 / 4 . Our analysis re- 
quires M 1 and for realistic lattice systems n < 5; con- 
sequently, the scaling in total atom number N = Mn is 
below the Heisenberg bound. Nevertheless, we find that 
measurements testing the predicted scaling behaviors for 
interaction-based metrologies should be possible in a lat- 
tice system using quadrature measurements. We note 
that measuring the quadratures of nonequilibrium lat- 
tice fields, our first goal, is less experimentally demand- 
ing than determining interaction strengths with optimal 
precision. Recently, super-Heisenberg scaling has been 
realized in another system [20] , 
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FIG. 1: (Color online) Schematic of quadrature operator in- 
terferometer for two-state atoms held in M optical lattice 
sites. Time evolution is left to right. The left-most box rep- 
resents initial state preparation. Solid and dashed lines corre- 
spond to atoms in lattice sites i and Bloch states with quasi- 
momenta k, respectively. Bold-red and thin-blue lines are 
atoms in spin states \a) and \b), respectively. Green lines in- 
dicate superposition states. Arrows specify spin states. Gray 
and red boxes correspond to unitary operations described in 
the text. Detectors (brown) at right measure population dif- 
ferences between spin states at quasi-momenta k. 



In the following, we first describe the interferometric 
technique and propose implementations. We next ana- 
lyze the special case of atoms held in a deep, cubic lattice 
with negligible tunneling, and obtain analytic predictions 
for the time evolution of both quadrature observables 
and their fluctuations. We emphasize, however, that the 
technique can provide time-resolved quadrature measure- 
ments for more complicated Hamiltonians and lattice ge- 
ometries (e.g., cubic, honeycomb, superlattices [2"5r[5T] ). 
including systems with tunneling, spin-orbit couplings 
using artificial gauge fields [S], and higher-band effects 
[32H33], and a new method for studying quench dynam- 
ics [35H37] . Related theoretical developments on collapsc- 
and-revival physics can be found in [10l I38H40] . 

Interferometry technique. Figure [T] shows a schematic 
of the quadrature interferometer. First, a superfluid of 
two-component bosonic atoms, with spin states \a) and 
| b) , are prepared in an optical lattice with all atoms in 
state \b). Interactions between &-state atoms [b atoms) 
are then slowly turned off with a Feshbach resonance, 
such that the initial interferometer state is a coherent 
state with k — quasi-momentum in the lowest Bloch 
band of an optical lattice potential. We therefore have 
|V>fc=o) = exp(7B^ =0 -7* J B fc=0 )|0), where |0) is the empty 
lattice, B k annihilates b atoms with quasi-momentum 
k, and the mean total atom number is N — |7| 2 . We 
can re-express this as the separable state \ipk=o) = 
11*=! exp(/36j - P*h)\0), where B k=Q = E^i/VM and 
bi annihilate b atoms in the lowest Wannier function of 
lattice site i. The mean atom number per lattice site is 
n — |/3| 2 = N/M. For notational simplicity, we suppress 
the band index and vectorial nature of i and k. As an 
aside, assuming an initial Fock state (Bl =0 ) N |0) /\/N\ 
(e.g., see [U]) leads to the same interferometer results 
found below with finite-size corrections of order A/~ 3 / 2 . 

The lattice depth is next increased to turn off tunneling 



without exciting atoms to higher bands, and the Wan- 
nier function can be approximated by the energetically- 
lowest local orbital in a site. At this point atoms nei- 
ther tunnel nor interact. A 7r-pulse between internal 
states a and b is then applied to half of the lattice 
sites, transforming the system state to nf=i 2 exp(/3a] — 
P*<h) ® UjLm/2+i ex P(/^j — P*bj)\0), where a { annihi- 
lates a state atoms (a atoms) in the lowest orbital of site 
i. The 7r-pulse is the first beam splitter of the interferom- 
eter, creating two arms with a and b atoms, respectively. 

After the 7r-pulse, atoms are held for a variable evolu- 
tion time t. The a atoms interact and undergo nonequi- 
librium dynamics, while use of the Feshbach resonance 
ensures that b atoms remain non-interacting and provide 
a well-characterized reference "beam" state. Spatial sep- 
aration between a and b atoms is required because it 
effectively turns off their mutual interactions. After this 
point, we can change the system dynamics to a "non- 
trivial" Hamiltonian, as long spatial separation of a and 
b atoms is maintained. For example, tunneling can be 
turned back on, the lattice geometry dynamically trans- 
formed, or external fields applied. 

The second beam splitter combines two operations. 
First, a Fourier transform on atom spatial modes is ap- 
plied, transforming site indices i into quasi-momentum 
indices k. This is achieved by releasing atoms from 
the lattice followed by (time-of-flight) expansion until 
all atoms have spatially overlapping wavcf unctions. Sec- 
ond, a 7r/2-pulse with tunable phase \ is applied to the 
internal states of all atoms, driving \a) —> (e lx ^ 2 \a) + 
e-^/ 2 \b))/V2 and \b) -> {-e tx ' 2 \a) + e~ lx / 2 \b)) / ^2. 

Finally, we measure the observable S k = A\A k — B k B k 
by detecting the difference between the number of a and 
b atoms in quasi-momentum k after the second beam 
splitter. Switching to the Heisenberg picture and using 
A k -> (e^/ 2 A k + e~^/ 2 B k )/V2 and B k -► {-e lx / 2 A k + 
e~' lx / 2 B k )/y/2 at the second beam splitter, we find that 
the observable is equivalent to S k —A\B k e~ tx +B\,A k e lx 
with expectation values taken with respect to wavefunc- 
tions just before the second beam splitter. 

Because the many-body wavefunction is separable in 
the internal atomic states before the 7r/2-pulse, it fol- 
lows that (S fc ) = {A k ){B\)e- lx + (A\){B k )e ix . The 
non-interacting reference atoms give (B k ) = f3g k with 
P = \ffie l6 and the complex-valued g k depends only on 
lattice geometry and distribution of b atoms. In partic- 
ular, g k= o — yM/2 for any lattice geometry and atom 
distribution. We therefore find 



(S k ) = 2Vn\g k \(X k ) 



(1) 



where X k is the quadrature observable with tunable 
phase Cfe = arg(e^ e+x ^fc). In contrast, regular collapsc- 
and-revival experiments |H HH 02], using atoms in 
only one spin state, measure the quasi-momentum distri- 
bution R k = A\A k . 
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Implementations. Realization of the interferometer re- 
quires two internal atomic states that have little or no 
two-body collisional loss, while allowing for tuning of 
b + b collisions to turn off interactions. We suggest 
the two energetically-lowest hyperfine states of bosonic 
alkali-metal atoms as only weak magnetic dipole-dipole 
and second-order spin-orbit interactions can cause inelas- 
tic loss. Candidates are 87 Rb near the narrow resonance 
at a magnetic field B of 100.7 mT 03] and 39 K near a 
much broader resonance at 40.2 mT [31], the latter re- 
quiring less stringent field control. Other examples can 
be found in [27J . Figure [2^, shows the scattering length 
a s of the 39 K resonance as a function of B and indicates 
two field values where one of the two hyperfine states has 
zero scattering length and, hence, does not interact. 

The initial state of the interferometer is prepared by 
starting with a (superfluid) Bose condensate of inter- 
acting 6-state atoms in a weak trap with level spacing 
huj vll typically harmonic with uj w /2tt — 1 Hz- 100 Hz; h 
is Planck's reduced constant. After adiabatically turning 
on the optical lattice, the system can be described by 
the single-band Bose-Hubbard Hamiltonian [45 4<3] with 
tunneling energy J and atom-atom interaction strength 
Ubb proportional to a s . In a simple cubic lattice, a su- 
perfluid ground state requires J 3> [7^/(5.83 x 6) at the 
end of the adiabatic ramp (the numerical factor follows 
from a 3D mean-field calculation). Typically, Ubb/h ~ 1 
kHz away from Feshbach resonances. Due to atom-atom 
interactions, at this stage the superfluid is only an ap- 
proximate coherent state. 

Figure (2}j shows the next step in initial state prepa- 
ration. Slowly turning off Ubb at fixed J using a Fesh- 
bach resonance over timescale Tint, the ground state ap- 
proaches the desired coherent state of non-interacting 
b atoms. Choosing r; n t 2ir/u> w minimizes excita- 
tion from zero quasi-momentum. After interactions are 
ramped-off, tunneling is suppressed by increasing lattice 
depth with time scale Tt un . Choosing T tun -C 2ir/(jjh g , 
where huj^ is the lattice band gap (typically w bg /27r = 10 
kHz-100 kHz), suppresses excitation to higher bands. 

The first beam splitter, a spatially selective 7r-pulse, 
can be implemented, for example, by illuminating only 
the top half of the lattice sites (see Fig. [2J:) with lasers 
driving optical Raman transitions via off-resonant tran- 
sitions to electronically-excited states. The interaction 
strength for a atoms is U aa . This approach is limited by 
the boundary over which the Raman laser intensities go 
to zero, which is no sharper than the laser wavelength. 
Uncertainty in the fraction of a versus b atoms leads to 
only small corrections of order M -3 / 2 . An alternative 
approach uses dynamically transformable lattices while 
simultaneously flipping selected spin states. For exam- 
ple, Refs. [Ml 02] used spin-dependent lattices to trans- 
form from a single- to double-well lattice where the spin 
state of atoms alternates with lattice site. In principle, 
alternating planes of a and b atoms can be created (see 
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FIG. 2: (Color Online) (a) Scattering lengths a s of 39 K col- 
lisions versus magnetic field B near a Feshbach resonance lo- 
cated at B w 40 mT. Curves are labeled by \fm) = |11) and 
1 10), the energetically lowest hyperfine channels of 39 K. Two 
operating points for the interferometer are shown by black- 
dotted lines. (Here lao = 0.0529 nm). (b) The interferometer 
initial state, a coherent state of non-interacting atoms, is cre- 
ated starting from a superfluid of interacting atoms by using a 
Feshbach resonance to slowly turn off the interaction strength 
Ubb, followed by faster turn off of the tunneling rate J. See 
text for time and energy scales. Panels (c) and (d) show two 
possible atom distributions after the site-specific 7r-pulse. In 
panel (c) atoms in state a (red spheres) and b (blue spheres) 
are located in the top and bottom half, respectively. In panel 
(d) planes of state a and b atoms alternate. 



Fig. [2ji) and our interferometric techniques used to study 
correlations in two-dimensional many-body physics. 

The second beam splitter involves release and expan- 
sion from the lattice followed by a global 7r/2 pulse on the 
internal state of all atoms. The latter operation requires 
no spatial selectivity and can be implemented with Ra- 
man transitions, or microwave and radio-frequency radi- 
ation. The phases 7r/2 and \ are controllable by choosing 
the polarization of the radiation relative to the applied 
magnetic field; the required techniques have been demon- 
strated in atomic clock and ion trap quantum computing 
implementations [49] . Requirements on expansion time 
are analyzed in |50j . and we estimate that ~ 20 ms is 
sufficient for operation of the interferometer. 

Finally, there are corrections due to interactions after 
atoms are released from the lattice. Following [30], we 
model the dominant effect by a time-varying U aa pro- 
portional to the inverse volume of the expanding local 
orbital. This leads to a small additional evolution time of 
order l/wb g , typically ~ 10 /zs, and independent of hold 
time t. Corrections due to interactions between atoms 
originating from different lattice sites are even smaller 
because of negligible overlap of their wavefunctions im- 
mediately after release. To minimize corrections to the a 
atom reference state and unwanted interactions between 
a and b atoms, the 7r/2 pulse should be delayed > l/wbg> 
but applied before significant phase shifts from magnetic 
field inhomogeneities arise. Effects of field gradients after 
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the 7r/2 pulse only lead to shifts in the relative location 
of the k = modes for a- and 6-atom clouds. 

Nonequilibrium dynamics. The factorization property 
leading to Eq. [T] holds even if tunneling and other mod- 
ified dynamics are induced for the interacting a atoms 
after the first beam splitter, making it possible to use 
the technique to measure the quadrature of nonequilib- 
rium fields for many types of Hamiltonians (e.g. systems 
with spin-orbit couplings, different lattice geometries, or 
applied fields) . Below, we analyze the interferometer and 
predict the quadrature evolution for the relatively simple 
evolution of the Bose-Hubbard model without tunneling, 
i.e. H MB = (U aa /2)J2t=i a i a laia l . Effective multibody 
interactions |51H53j are omitted for simplicity. 



We focus on the normalized k — quasi-momentum 
observable s (t) = (Sk=o(t))/N, which is proportional 
to the order parameter. (For other Hamiltonians k ^ 
quasi-momenta may be interesting.) Defining single-site 
quadrature operators Xi = (e~^°ai + e"*°al)/2, we find 
s o(t) = (xi(i)) I (2y/n), after expanding the momentum 
operator in terms of site operators and realizing that the 
time evolution is the same for all sites. For Hmb we 
derive (ac<(*)) = \/Hexp([cos(<^) — l]n) cos(n sin(</>) + Co), 
with phase (ft = U aa t/h. Figures [3^, and b show exam- 
ples of the periodic evolution of (xi(t))/\/% for several 
values of Co an d ^- The interferometer makes possible 
direct experimental verification of these quadrature dy- 
namics. The time evolution ofTZk—o/N = \(a,i)\ 2 /n (as- 
suming large M), determined by regular collapse-and- 
revival measurements [U rT2TTH] , is also shown. 

The fluctuations in quadratures observable for 
nonequilibrium Bose gases in deep lattices are also in- 
teresting. Figures |3J: and d shows the analytic time- 
dependent (single-site) uncertainty a X i — ((acj — (xj)) 2 ) 1 ' 2 
as a function of time. For a coherent state for any n and 
Co, Cxi = 1/2, corresponding to the SQL and shown by 
horizontal dashed lines; a X i < 1/2 indicates quadrature 
squeezing. For most times and phases Co, the uncertainty 
is a xi > 1/2; in fact, a "plateau" develops for large n 
where a X i approaches \/2n + 1/2. For times near integer 
multiples of U aa t/h and values of Co near zero, however, 
a X i < 1/2, and the fluctuations in (x^) are reduced below 
the SQL. This is seen in Fig. |5J: for Co = 0, where we 
find that the minimum value of a xi approaches a con- 
stant ss 0.29 for large n. In contrast, Fig. [3}i shows that 
a x i for Co — 7r /2 is always antisqueezed. 

Interaction-based metrology. In addition to directly 
measuring field quadratures, the interferometry can be 
used to measure interaction strengths, although achiev- 
ing optimal performance is experimentally more demand- 
ing. The period h/U aa , and hence U aa , can be obtained 
by determining times t\ and t 2 > t\ at which so(t) and 
its time derivative are the same. Figure [3^, shows an ex- 
ample where U aa = /i/(<2 — £i). Error propagation gives 
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FIG. 3: (Color online) Panels (a) and (b) show the normal- 
ized single-site quadrature {xi)/^/n (solid lines) as a function 
of evolution time in units of h/U aa for various quadrature 
phases and mean atom number n. The double-sided arrow 
indicates a possible time interval for measuring U aa . The den- 
sity operator \{ai)\ 2 /n (dashed line) is only shown for n = 4. 
Panels (c) and (d) show the quadrature uncertainties a X i as 
a function of evolution time and n, for = and Co = tt/2 , 
respectively. The dash-dotted line, labeled SQL, indicates 
&xi = 1/2 for a coherent state. 
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with expectation values evaluated at phase 4>i = U aa t\ / H, 
assuming equal contributions at t\ and ti added in 
quadrature. The fractional uncertainty gives SQL scaling 
oc 1/ y/M/2 in number of a-atom lattice sites, as expected 
for a site-separable wavefunction, and we can view the in- 
terferometer as M/2 independent probes of n interacting 
atoms. The second term under the square root is due 
to the uncertainty in total atom number N. The un- 
certainty also scales as 1/m in the number of oscillation 
periods m; for brevity we set m = 1. 

Minimizing 6U aa /U a a based on measurements of Xi(t) 
involves a trade-off between identifying optimal measure- 
ment times ti, ti (or phases <j>\,4>2) that maximize the 
slope d{xi)/dcp while minimizing a X i. We find that the 
smallest (optimal) fractional uncertainty 



(SU aa /U aa ) queid = (l/2 7 r)M- 1 / 2 n- 3 / 2 , 



(3) 



is obtained from the quadrature observable with Co = 
7r/2, and measurements made when (x^ — and the 
slope d(xi) /defi is maximum. This occurs at phases 4>i = 
and 4>2 — 27r, corresponding to times t\ — and t 2 — 
2irh/U aa . Because (x^ = and a X i — 1/2 at these times, 
the result is independent of uncertainty in N. For other 
Co, the numerical prefactor in Eq. ^ is larger, but the 
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scaling in n is the same. As noted in Implementations, 
the small correction from interactions during expansion 
can be included in the determination of tip- 

At optimal measurement times the state is unsqueezed, 
and the super-Heisenberg scaling n~ 3 / 2 derives from the 
nonlinear dependence of d(xi) /dip on n. Our quadra- 
ture interferometer yields the best possible scaling be- 
havior predicted for interaction-based metrology schemes 
based on two-body interactions and non-entangled states 
[16j . We have also analyzed the minimal uncertainty 
possible from conventional collapse and revival mea- 
surements of (1Zk=o(t)) / N . We obtain the analytic 
estimate (SU aa /U aa ) coav = (2/3) 3 / i n- 1 M- 1 / 2 n- 3 / i , 
and have confirmed it is in good agreement with nu- 
merical evaluation. The scaling is significantly less ad- 
vantageous than what is possible using the quadrature 
interferometry. 

We can compare our predictions with conventional 
collapse-and-revival experiments. Will et al. [12] ob- 
tained SU aa /U aa ps 2 x 1CT 2 with N = 2 x 10 5 atoms 
and n = 2.5. For these parameters (SU a a/U aa ) conv — 
3 x 10~ 4 . The larger experimental uncertainty is likely 
due to fluctuations in the total atom number and lattice- 
laser intensities. For the quadrature interferometer 
(5U aa /U aa ) quad = 1.2 x 1CT 4 is, in principle, possible, 
giving nearly a factor of three improvement even with 
h = 2.5. As noted above, we also predict less sensitivity 
to fluctuations in N using the quadrature technique. 

Outlook. Our analysis suggests that even for small (re- 
alistic) n it should be possible to test predicted scaling 
behaviors for interaction-based metrologies in an opti- 
cal lattice system. There are many interesting directions 
for future research. For example, are there interferome- 
ter states that give n~ 2 scaling per lattice site, as sug- 
gested by the analysis in |16| ? Are there implementa- 
tions where interesting information lies in k ^ quasi- 
momenta states? Can we modify the interferometer to 
measure observables such as (a,Oj), and use entangle- 
ment between atoms in different sites to measure nonlo- 
cal field correlations? There are also practical challenges 
facing implementations, including the effects of inhomo- 
geneities, lattice (laser) fluctuations, and limits on both 
atom number per site and total atom number. Even with 
non-optimal measurements, however, the interferometer 
technique developed here should expand the toolkit for 
studying nonequilibrium dynamics in optical lattices. 
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